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Abstract 

The lattice spin model with Q -com P onen t discrete spin variables restricted to have 
orientations orthogonal to the faces of Q-dimensional hypercube is considered on the Bethe 
lattice, the recursive graph which contains no cycles. The partition function of the model 
with dipole-dipole and quadrupole-quadrupole interaction for arbitrary planar graph is 
presented in terms of double graph expansions. The latter is calculated exactly in case of 
trees. The system of two recurrent relations which allows to calculate all thermodynamic 
characteristics of the model is obtained. The correspondence between thermodynamic 
phases and different types of fixed points of the RR is established. Using the technique of 
simple iterations the plots of the zero field magnetization and quadrupolar moment are 
obtained. Analyzing the regions of stability of different types of fixed points of the system 
of recurrent relations the phase diagrams of the model are plotted. For Q < 2 the phase 
diagram of the model is found to have three tricritical points, whereas for Q > 2 there 
are one triple and one tricritical points. 

PACS number (s): 05.50.+q, 05.70.Fh 







1 Introduction 



Cubic symmetry plays a significant role in many types of phase transitions and critical phenom- 
ena. Still in 70-s various models possessing cubic symmetry have been introduced to determine 
the nature of the displacive phase transitions in perovskites [1, 2]. Magnetic phenomena in 
cubic crystals are also affected by the lattice structure. For instance, in the crystalline solids 
with cubic-symmetric lattices (Fe, Ni, etc.) it leads to the modification of the magnetic ex- 
change interaction giving rise to additional contributions to the conventional 0(N)-symmetric 
Heisenberg Hamiltonian. The simplest contribution of the underlying cubic symmetry is the 
single-ion anisotropy of the form J2i s t ■ During the last three decades different aspects of this 
issue have been the subject of various investigations. 

In the framework of the field-theoretical approach to critical phenomena modelling of effects 
of cubic symmetry is usually performed in terms of continuous-spin Landau-Ginzburg effective 
Hamiltonian with cubic anisotropy, i.e. the 4 - theory with an additional cubic term which 
breaks explicitly the O(N) invariance to a residual discrete cubic symmetry [3] - [7]: 



where (f>(x) = (</>i,<f>2, ■ ■ ■ ,</>n) is a continuous N-component local vector order parameter, 
4 = (0 2 ) 2 and P — 1/kT is inverse temperature. The spin orientations, orthogonal to faces of 
an iV-dimensional hypercube will be favorable for v < 0, whereas positive values of v favor the 
orientations toward the corners. 

This model received numerous important applications, among which are such as the oxygen 
ordering in YBaCusOe+x [8], one of the most studied high-T c superconductor; the buckling 
instability of confined colloid crystal layer [9]; the micellar binary solution of water and am- 
phiphile [10] etc. The two-dimensional case has been recently intensively studied within the 
framework of renormalization group technique in the space of fixed dimensionality up to five- 
loop approximation [5]- [7]. The model was found to have four fixed points: the Gaussian one, 
the Ising one with N decoupled components, the 0(N)-symmetric and the cubic fixed point. 
Analogous calculations were made for D = 3, as well revealing some peculiar type of cubic fixed 
point for N > 2 corresponding to specific anisotropic mode of the critical behavior [11]-[15] . 
The surface critical behavior of the model was also studied [16]. It has been recently shown 
that the behavior of the spherical many-spin magnetic nanoparticle with surface anisotropy 
may be modelled by an effective single macro-spin with cubic anisotropy terms in the effective 
energy [17]. 

Another class of lattice cubic-symmetric spin models was introduced by Kim, Levy and 
Uffer in Ref. [18] to explain the tricritical-like behavior of cubic rare-earth compounds, partic- 
ularly holmium antimonide, HoSb [18]-[23]. Projecting the pair exchange interaction onto the 
sixfold degenerated ground-state manifold of Ho 3+ ion they arrived at the following effective 
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Hamiltonian: 

-PH = '£j ij S i S j , (2) 

(id) 

where the classical spin variables S, are the unit vectors restricted to have orientations, orthog- 
onal to the faces of the cube and the sum is going over all the pairs of nearest-neighbor sites. 
Generalization for Q component spin is obvious: each spin can assume 2Q orientations: 

S t e {(±1,0,...,0),(0,±1,...,0)...(0,...,0,±1)}. (3) 

This model is known as the Face-cubic model. It is connected to the continuous cubic model 
of Eq. (1) via the limit of strong anisotropy (\v | ^> |6|). In a similar way one can also 
consider quadrupolar pair interaction terms and obtain the following Hamiltonian provided all 
interactions are homogenous [20, 21]: 

-PHfc = J E( s * s i) + K E(S^-) 2 . (4) 

It is easy to see that Hamiltonian (4) can be represented in terms of two sets of discrete variables: 
a Potts-like determining which component of Sj is non-zero and an Ising-like, corresponding 
to the sign of the component. Indeed, 

SjSj = O'iO'jSaiaji (5) 

where <j{ — ±1 and = 1, 2, ...Q. Therefore we have 

-PH F C = J E a i a + K E S <*i,a r (6) 
ihj) (hj) 

Formally we can enlarge the Hamiltonian (4) up to the interaction terms of higher power of the 
(S.S,): 

-mf c = E E USiSiT- (7) 

(ij) n=0 

For arbitrary finite L this expression leads to the same Hamiltonian that of Eq. (6) with 
J = J2k=o ^2k+i an d K = J2k=o J^k- At K = the Hamiltonian (7) reduced to the 2Q - 
state Potts model; J = corresponds to two decoupled Q - state Potts model, and at Q = 2 
one obtains the Ashkin-Teller model. Variational renormalizat ion-group study of the pure 
and diluted Q - component face-cubic model in two dimensions has revealed existence of four 
competing possible types of critical behavior corresponding to Q - state Potts model, 2Q - 
state Potts model, Ising model and special "cubic" fixed point [22, 23]. It was also found that 
at Q < Q c = 2 the transitions are continuous and critical behavior of the discrete face-cubic 
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model belongs to the 0(N)-model universality class. For Q = 2 the Ashkin-Teller-like behavior 
occurs and for Q > 2 transitions are of first order. There are also early results obtained within 
the mean-field (MF) theory [18] and using the Bethte-Peirels (BP) approximations and high- 
temperature series [19] for the case when only dipolar pair interactions are included (K = 0). 
In the BP approximation the critical value of spin component Q c , above which transitions are 
of the first order was found to be given by 



where q is the coordination number of the lattice. The limit q — > oo corresponds to the MF 
solution and gives Q c = 3. In Ref. [19] the high-temperature series for Q - component face- 
cubic model on three dimensional fee lattice were constructed up to 5-th order, from which 
authors obtained Q c = 2.35 ± 0.2, whereas Eq. (8) gives us Q c ~ 2.8 Moreover, although the 
MF solution of cubic model predicted the tricritical like behavior [18], the BP approximation 
showed that it is not the case [19]. Only inclusion of single-ion-anisotropy terms, quadrupolar 
pair interactions and crystal fields may drive the system tricritical [20]. 

The paper is organized as follows. In Section II we give the introduction to recursive lattice 
constructions and their role in statistical mechanics. Then we discuss the peculiarities of the 
Cayley tree, Bethe lattice and connections between these two objects. Section III is devoted 
to the graph expansions for the partition function of the model under consideration defined 
on the arbitrary planar graph. We presented the partition function in terms of double power 
series which in case of the Cayley tree are summed up exactly, giving the same expression as 
in case of one-dimensional linear chain. We also discuss arguments of Eggarter concerning the 
mechanisms of phase transitions on the Cayley tree and Bethe lattice. In the next Section 
IV we derive the system of recurrent relations for the FCg-model on the Bethe lattice refor- 
mulating the statistical-mechanical problem in terms of the theory of dynamical systems [25]. 
We introduce the order parameters and give their expressions in terms of recursive scheme. In 
Section V we identify different types of the fixed points of the system of recurrent relation with 
different physical phases and present the plots of order parameters temperature dependencies. 
The phase diagrams for different values of spin component number Q are also presented. Brief 
summary is given in Section VI. 



2 Recursive methods, Cayley tree and Bethe lattice 



Among the vast variety of statistical mechanics lattice models with strong local interactions 
only very limited amount allows exact solutions. These exact solutions are known only for 
low-dimensional systems, more precisely for d = 1 and d = 2; most of exact solutions for two- 
dimensional systems are known only in the absence of the external field and/or for the special 
choice of model parameters [24] . Well known conventional approximate methods like MF theory 
and BP approximation in general can provide only more or less qualitatively satisfactory picture 
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Figure 1: The Cayley tree with coordination number q = 3 and 3 shells. 



and in some cases they just fail. That is why the quest for the alternative approaches, which 
can provide more reliable results for the thermodynamics of lattice models is very important. 

The recursive lattices are of twofold interest. On the one hand, the models, defined on a 
recursive lattices can be considered as an independent original area of research in which the 
powerful methods of dynamical systems theory and fractal geometry are successfully exploited 
for determining the thermodynamical properties of the statistical mechanics models. On the 
other hand, recursive lattices provide a specific kind of approximate treatment of many-particle 
systems physics. In the heart of these approximations not the simplification of the character 
or/and strength of the interaction between the system elements, but the modification of the 
topology of the underlying lattice lies. This modification in most cases consist in the replace- 
ment of the regular periodic "physical" lattice by a recursive one, constructed by the certain 
algorithm and possessing the self-similarity, maintaining all interaction unchanged. For the 
large class of lattice models with commuting variables (" classical models" ) this approach leads 
to the exact solution of the statical mechanics problem in terms of the theory of dynamical 
systems or, more precisely, in terms of discrete maps. The solution formally is rather similar to 
that which corresponds to conventional BP approximation but, in contrast to the latter, it is 
the exact solution for the lattice endowed with the recursive structure provided the boundary 
sites are properly taken into consideration. It was argued that in some cases, particularly in the 
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models where multi-site interactions are presented [31]-[33] the recursive lattice approximation 
gives more reliable results than conventional BP [34]. 

The simplest example of recursive lattices is a Cayley tree, the self-similar graph, which 
contains no cycles. In order to construct a Cayley tree one should start with a central site O. 
At the first step this site should be connected by links with q others, which constitute the first 
shell or the first generation of the Cayley tree. The second shell of the Cayley tree is constructed 
by repeating this procedure for all sites of the first shell: each of the q sites belonging to the 
first shell connects to the q — 1 new sites. Thus, the second shell contains N 2 = q(q — 1) sites. 
Accomplishing this construction for n steps one arrives at the recursive connected graph, which 
contains no cycles and is called the Cayley tree with coordination number q and n generations 
(shells) (Fig. 1). As the number of sites in k-th shell is q(q — l) k ^ 1 the total number of sites 
then 

Af. = l + t i(l-lf- 1 = q(q ~' r ~ 2 - (») 

k=l 1 Z 

Thus, each site of the Cayley tree has coordination number equal to q except the sites on the 
last shell which have only one neighbor. The peculiar property of the Caylet tree is the large 
amount of boundary sites which even in the thermodynamical limit n — > oo comprise a finite 
fraction of the total number of sites 

Hm ,(«-!)- = Mm fa (j ; = 

n^oo N n n^oo q (g _ l) n _ 2 q-\ V ' 

This feature causes anomalous properties of the statistical mechanics systems on the full Cayley 
tree. For instance, in Ising model there is no zero field magnetization, whereas the derivations 
of free energy with respect to the external field exhibit singular behavior [26] - [29]. 

The Bethe lattice is an object intimately linked to the Cayley tree but, in contrast to the 
latter, it is devoid of the complications which arise from the boundary sites. When one deals 
with the Bethe lattice the general structure and topology of the Cayley tree are preserved 
whereas only contributions from the bulk sites lying deep inside the Cayley tree are taken into 
consideration. In this regard the Bethe lattice is the interior of the Cayley tree [24, 30]. Dealing 
with the Bethe lattice we suppose the underlying Cayley tree to be large enough to achieve 
thermodynamical limit and consider only the sites situated far away from a boundary. One can 
regard these sites as a uniform lattice with coordination number q. 

An undeniable advantage of the Bethe lattice is the possibility of exact solutions for many 
types of statistical mechanics problems. A large amount of problems of the theory of magnetism, 
macromolecule physics, lattice gauge theory, self-organized criticality, dynamical mean-field 
theory (DMFT) and general questions of statistical mechanics like zeroes of partition functions 
have been successfully considered on Bethe or Bethe-like lattices, revealing many interesting 
exact results and deep connections between the theory of dynamical systems and statistical 
mechanics [26] - [49]. 
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3 Face-cubic model on planar graphs 



It is possible to calculate exactly the partition function of the FCQ-model on Cayley tree in 
case of the absence of external field. Let us consider the partition function of the model without 
external field on arbitrary planar graph G 

^ = EE IT exp{(Ja t a J + K)5 ai , aj }. (11) 
M M (*J)eG 

As in the case of Potts model [50] one can represent it in the following way: 

z g = EE II (i + tWc^.), (12) 

M {a} (iJ)eG 

where U CiGj = e J(Tl<T: > +K — 1. Then summing out over the all we obtain that 

% = EE^ (ff) II u*w (is) 

M G'CG (i,j)eG> 

where the second sum is going over all spanning graphs G' of the underlying planar graph G 
and Rq(G') is the 0-th Betti's number of G 1 which coincides with the number of connected 
components of G' and the product is over all links belonging to G' . Using the identity 
(TjCTj = 25 ai<7j — 1 one can represent Eq. (13) in terms of double power series associated with 
the underlying lattice provided latter is a planar graph: 

z G = E Q %(G V (G,) E n (i + f<W,)> 

G'CG {a} (i,j)eG> 

u = e K ~ J -l, 
e 2J - 1 

v = T^k' ( 14 ) 

where e(G') is the edges number in the graph G' and all Oi assume values or 1. So, for each 
spanning subgraph G' of planar graph G we obtain a partition function of 2-state Potts model: 



i G = Q Ro{G,) u< G,) z?? tts {v) 

G'CG 

= e Q Ro(G '^ e(G,) E 2^ G V G "\ (is) 



G'CG G"CG' 



where the second sum is going over all subgraph of G' Thus, we have succeeded in representing 
the partition function of the FCq - model given by Hamiltonian (6) in terms of double power 
series associated with the underlying lattice provided the latter is a planar graph. As in case 
of ordinary spin Potts model one can make 1/Q-expansions for FCQ-model with the aid of Eq. 
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(15). For instance, for the two-dimensional square lattice one obtains 

Z = (2gni + 2iV^ + l) + Ar(2iV-l)(^) 2 (| + l) 2 + C 2 3 7V (^) 3 (^ + l) 3 (16) 
+ {Ct N - N)(^n V - + l) 4 + N^(( V - + l) 4 + C-f) + (C% N - N(2N - 4))(^) 5 (| + l) 5 

+ 2iV^((| + l)^ + 2((|r + (^) + 3(|) 6 ) + ...}. 

Now we can use the well known relations between the topological invariants of graphs [51] to 
rewrite Eq. (15) in another way suitable for our further purposes. Applying the Euler theorem 
to planar graph G we get 

R (G)-R l (G) = n-e(G), (17) 

where the first Betti's number R±(G) coincides with the number of independent cycles of G 
and n is the number of sites. Therefore 

Z G = (2Q) n E Q^(G')^y(G>) £ 2 R ^ G "\-)< G "\ (18) 

G'CG Q G"CG' ^ 

3.1 The case of graphs without cycles: trees and forests 

Partition function for FCg-model in the form of Eq. (18) can easily be calculated in case of 
the so-called forests, the planar graphs, which contain no cycles. Connected part of forest is 
called tree. Let us suppose that we have a tree T with n sites and L n edges. Then Eq.(18) 
takes the form 

Zt = {2QT E O e{G ' ] E & {G " ] 

G'CT y G"CG' 

= (2QrE^(^) fc E^(^ (19) 

k=0 l=o z 

where C l k = ^A, are the binomial coefficients and the sum is going over all subgraphs with 
given volume k. With the aid of Newton's binom we obtain 

M «"(S(f +i ) +i f ■ (20) 

If T is the Cayley tree then L n = n — 1 and we will have 

Zcayley = 2Q (u (v + 2) + 2Q)"- 1 . (21) 
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Thus, the free energy per site for FCQ-model on Cayley tree is 



/ = -k B T lim 



k B T log 2 - k B T log (e K cosh J + Q - l) . 




(22) 



It is noteworthy that this result is in full agreement with that obtained by Aharony for 
one-dimensional model by transfer-matrix technique [21]. Indeed, a one-dimensional chain 
can be regarded as a "tree" in the sense mentioned above. Thus, the free energy of the FCq 
model on a Cayley tree in a thermodynamic limit is continuous for all T. Formally it coincides 
with that for one-dimensional systems where the continuity is really the case. But, as is known, 
the BP approximations for the lattice model with an arbitrary number of component solves 
exactly the problem on a Cayley tree [40] . Therefore the model on a Cayley tree must possess 
a finite critical temperature as predicted by BP approximation [19]. The origin of this peculiar 
properties of the Cayley tree was revealed by Eggarter [26] on the example of Ising model. He 
argued that the equivalence of all sites in the thermodynamic limit, which is one of the key 
points of BP approximation, breaks down for the sites of the Cayley tree which are situated 
close to the surface. As appears from Eq. (10) these sites comprise rather large fraction of the 
sites of the Cayley tree and, thus, they determine the behavior of all thermodynamic quantities 
to a great extent. 

However, even though no phase transitions in the thermodynamic limit occur at the value 
of critical temperature T c predicted by the BP method the value of the order parameter (mag- 
netization etc.) in the interior of the lattice, i. e. for the region far from the surface, which 
is generally called "Bethe lattice", will undergo jumps from zero to its BP value when the 
temperature is below T c . In order to illustrate it we will use the technique of dynamical system 
theory, which becomes a powerful tool for investigating various physical problems on recursive 
lattices. 

4 Recursive method for Face Cubic model on Bethe lat- 
tice. 

Many statistical systems defined on the recursive lattices are famous for the possibility of 
exact solution in terms of dynamical system theory. In the heart of these exact solutions 
the self-similarity of recursive lattices lies. For instance, if we cut apart the Cayley tree at 
the central site it will give q identical branches each of which is the same Cayley tree with 
the number of generation decreased by 1. Using this fact one can establish the connection 
between the partition function of a model defined on the Cayley tree containing n shells with 
the partition function of the same model defined on the Cayley tree containing n-1 shells. 
Therefore, the thermodynamical problem is reformulated in terms of discrete maps, given by 
recurrent relations. Let us introduce a symmetry breaking field into Hamiltonian (6). The 
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field can be considered as a uniform magnetic field pointing along the first coordinate axis 
H = (H, 0, ...0). The corresponding interaction Hamiltonian is of the conventional Zeeman 
type: 

H Z = -UJ2S 1 , (23) 

i 

which leads to the following form of the FCg-model Hamiltonian in the field: 

-PH F C = J Wjtianaj + K fiauaj + h £ 0^,1, (24) 

{hi) * 

where h = f3H . Hereafter we pass from the Cayley tree to the Bethe lattice having in mind 
one of the ways mentioned in Section II. Thus, we should no longer care about boundary sites 
and boundary conditions. According to that one can represent the partition function of the 
FCg-model on the Bethe lattice in the following form: 

Z= £ e^M<?nK«o)] 9 , (25) 

where a and a are the variables of the spin in the central site and g n (a , a ) refers to a 
partition function of individual branch: 



g n (a , a ) = ex P ( Ja oVi + K) 5 aom + £ (Ja^j + K) 8 auaj + £ (Ji5 ai:1 . (26) 

Each branch, in its turn, can be cut at the site which was previously connected to the central 
one. This will give us q-1 identical branches being the Bethe lattices with n-1 generations. Thus, 
the connection between g n and g n _i is 



g n (a , a ) = £ ex P + K) 5 aQiCei + ha^^) \g n ^{a u «i)] 9 \ (27) 

(ci,q;i) 

Here we obtain the system of 2Q recursion relations but, in fact, only three of them are 
independent as all g n (a, a) for a ^ 1 are identical. So, from (27) we have 



g n (+, 1) = e J+K+h [g n ^(+, l)]*" 1 + e - J+K ' h [g n ^(-, l)]^ 1 + 2(Q - l)[g n -i(±, (28) 

9n(- 1) = e- J+K+h [g n ^(+, l)]*- 1 + e J+K - h [g n ^(-, l)]^ 1 + 2(Q - l)k-i(±, *)]«-\ 

g n (±, *) = e h [g n ^(+, l)]^ 1 + e^g^-, l)]^ 1 + (e J+K + e- J+K + 2(Q - 2))[< ? „_ 1 (±, *)]«" \ 

where g n (±, *) stand for any 2Q-2 partition functions corresponding to the individual branch 
with Si whose direction is non collinear with the first coordinate axes. Introducing the variables 

x n = g n (+, l)/g n (±,*), 

y n = g n (-l)/9n(±,*), (29) 
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we obtain the system of two recurrent relations 



%n fl\%n— 1) Un— l)j 

Vn = f2(Xn-l,y n -l), (30) 



with 



Pi(x,y) a^- 1 + 6/x- 1 ^- 1 + 2(Q - 1) 

i?(x,y) ^-i + ^V- 1 + a + 6 + 2(g-2)' 1 ' 

f( s P 2 (x,y) 6/i^ 1 + a^V 1 + 2(g-l) 

mz, y) - 



i?(x,y) /i^- 1 + / u-V" 1 + a + & + 2(Q-2)' 
where the following notations are adopted 

a = exp(J + K), (32) 
b = exp(-J + K), 
\x = exp(h). 

In this approach statistical averages of all physical quantities can be expressed in terms of x 
and y variables defined by Eq. (30). For instance, when the total number of spins is n the 
magnetization along the first coordination axis which is the thermal average of the following 
form: 

in i n 

"» = -E( s f ) ) = -E(*. 1 )). ( 33 ) 

in terms of x and y are expressed as 

fix q - //-y 

"1 = : 7TZ T, 34 

lixi + + 2(Q - 1) V ; 

provided all sites of the lattice are equivalent. It is easy to see that this quantity plays the role 
of an order parameter of Ising universality class. Another order parameter belonging to the 
Q-state Potts universality class is 

1 n 

p = -EW«i» 1 )>» (35) 
n i=i 

which also can be regarded as a "quadrupolar" moment (Si ). Thus 

= ^ + ^~ lyq (16) 

P ^ + /x-V + 2(Q-l)' 1 ' 
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The order parameters m and p define the three possible phases of the FCq model in case 
of the ferromagnetic couplings, i.e. J > 0, K > 0, provided external magnetic field is vanished: 



(a) disordered (paramagnetic) phase: m — 0, p — 1/Q 

(b) ferromagneticaly ordered phase: m ^ , p ^ 1/Q 

(c) partially ordered (quadrupolar) phase: m — , p ^ 1/Q 

5 Investigations of phase transitions in terms of dynam- 
ical systems theory 

In order to obtain physical results one must implement an iterative procedure for the RR (31). 
Namely, starting from the random initial conditions (xq, yo) one uses the simple iterations 
scheme and examines the behavior of physical quantities after a large number of iterations [35] . 
In a simplest case the iterative sequence {x n ,y n } converges to a fixed point (x*,y*), which is 
defined by 



This situation is inherent to the ferromagnetic case when both J and K are positive. As can 
easily be seen from Eqs. (33) - (36), when the magnetization and quadrupolar moment of the 
system take values m and p respectively, then the corresponding fixed point of the RR (31) is 
the following, provided h — 0: 



Thus, according to the properties of different phases of the ferromagnetic FCq model one can 
establish the connection between them and the classification of the possible types of fixed points: 

(a) disordered (paramagnetic) phase: x* = y* = 1 

(b) ferromagneticaly ordered phase: x* ^ y* ^ 1 

(c) partially ordered (quadrupolar) phase: x* = y* ^ 1 

The condition x* = y* leads to g n {+, 1) — 9n(~, 1)) which means that the probability of spin 
"up" is equal to the probability of spin "down" which is really the case in paramagnetic phase. 



x* = fi(x*,y*) 
y* = f2(x*,y*) 



(37) 




(38) 
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Figure 2: The temperature behavior of the order parameters for Q = 2 and q = 3. The transitions from 
disordered to ferromagnetic phase is of second order for all range of K .(a) Zero field magnetization for 
K/J = 0; 1 and 2.5; (b) Quadrupolar moment in zero field for K/J = 0; 1; 1.4 and 2. For K/J = 1 one 
can see the first order transition between disordered and ferromagnetic phases, whereas at K/J = 1.4 
one can see two subsequent phase transitions, the continuous transition from disordered to quadrupolar 
phase and very close to that the first order transition to the ferromagnetic phase. With the further 
increase of K/J the second transition becomes of the second order. 



Applying the simple iterative scheme to Eqs. (30), (34) and (36) one can obtain the plots of 
magnetization processes (m vs. H at fixed values of T) as well as the zero field magnetization 
and quadrupolar moment of the system. 

In Fig. 2 one can see the temperature dependencies of the order parameters of the FCq- 
model on Bethe lattice with coordination number q = 3 and Q = 2. As is obviously seen from 
Fig. 2(a) in this case the zero field magnetization is always continuous, whereas the behavior 
of p is quite different. Depending on the value of ratio Kj J the p(T) curve could be continuous 
or can include one first order transition point. In Fig. 2(b) one can see that at small values 
of K/ J only the one phase transition occurs from disordered phase to ferromagnetic phase 
omitting the partially ordered phase. This transition is of the second order until K/ J reaches 
some intermediate value above which one can see a discontinuity in the order parameter p 
at critical temperature corresponding to the transition between disordered and ferromagnetic 
phases. At large values of K/ J the system undergoes successively two phase transitions at 
temperatures T q {K) and Tf(K), the larger one corresponds to the transition between disordered 
and quandrupolar phases, the lower one - to the transitions from quandrupolar to ferromagnetic. 
Within this region of the values of K one can see the second order transition from disordered 
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j : 
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1.25 - 
1 - 



0.5 1 1.5 2 2.5 

K 

J 

Figure 3: The phase diagram of the model for Q = 2 and q — 3. Solid line corresponds to the 
second order transitions, dotted - to the first order transitions. For small values of the ratio 
Kj J one can see only single phase transition line between disordered hight-temperature phase 
and completely ordered ferromagnetic phase, whereas for K/J > 1.233 there is the partially 
ordered "quadrupolar" phase between them. The tricritical points are K/J 0.32; 1.233 and 
1.6. 

phase into "quadrupolar" phase and the further first order transition to the ferromagnetic 
phase. However, beginning with some value of K the transition between quadrupolar and 
ferromagnetic phases becomes continuous. 

In order to complete the picture of the phase structure of the model under consideration at 
Q = 2 we plot a phase diagram by separating the regions of the fixed points of different kinds 
of the RR from Eq. (31) in the (K/ J, Tj J)-plane. One can see in Fig. 3 that in case of Q — 2, 
q = 3 at K = the phase transition from disordered to ferromagnetic phase is of the second 
order. This feature maintains up to K/J m 0.32, where one can notice the tricritical point 
separating the region of continuous transitions from the region of the first order transitions. 
The line of the second order phase transitions between disordered and quadrupolar phases 
merges the line of transitions between disordered and ferromagnetic phase at K/J 1.233, 
the latter, in its turn, again becomes of the second order at K/J as 1.6 and for large values of 
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Figure 4: The temperature behavior of the order parameters for Q = 3 and q = 3. Here magnetization 
of the system undergoes a jump for small values of K/J and becomes continuous beginning with 
K/J ~ 1.84. Another order parameter always remains discontinuous at the transition points. In the 
plots p(T) for K/J = 1.7 one can see that the transitions between " quadrupolar" and ferromagnetic 
phases is of the first order, whereas with the further increase of K/J in becomes continuous. 



Kj J it becomes parallel to the Kj J axis which means that the transition temperature between 
quadrupolar and ferromagnetic phases is unaffected by the value of K beginning with Kj J pa 4. 

The situation is quite different for Q > 2. In Fig. 4 the plots of zero field magnetization 
and quadrupolar moment for Q = 3 and q = 3 are presented. In this case both m and p are 
discontinuous at the critical temperatures. Moreover, this feature is preserved for p for all values 
of K/J, whereas m becomes continuous after K/J pa 1.84. the corresponding phase diagram 
is presented in Fig. 5. Here one can see a triple point at K/J pa 1.67 where the merging 
of the three lines of the first order transitions takes place. There is also a tricritical point at 
K/J pa 1.84 which separates the region of the first and the second order phase transitions on 
the line between quadrupolar and ferromagnetic phases. The line, as in case of Q = 2, goes 
to plateau beginning with K/J pa 5, so the transition temperature again does not depend on 
K for K > 5J. In general, for Q > 3 the pase diagrams of the model under consideration 
have the same topology as in case of Q = 3. The main feature is the decrease of the distance 
between triple and tricritical points with the increase of the Q. So, for Q = 30 one can obtain 
that this distance is of 10~ 2 order with K/J. Apparently, they have never shrunk to a single 
point at finite values of Q, but it could be the case when Q — > oo. The value of critical spin 
component number Q c at which the phase transitions become of the first order at K = in the 



14 




Figure 5: The phase diagram of the model for Q = 3 and q = 3. Solid line corresponds to the 
second order transitions, dotted - to the first order transitions. The inset shows the merge point 
of the three first order lines (triple point): disordered-ferromagnetic, disordered-" quadrupolar" 
and "quadrupolar"-ferromagnetic. In this point at K/J as 1.67 these three phase are in the 
equilibrium . The tricritical point on the line between " quadrupolar" and ferromagnetic phases 
is situated rather close to the triple point. The distance between them tends to zero with the 
increase of Q. 



15 



Face-cubic model on Bethe lattice is exactly equal to 3, which coincides with the mean-field 
results at q — > oo, but, in contrast to that, on Bethe lattice this value is unaffected by the 
value of coordination number q. Obviously, this is the direct consequence of the fact, that the 
Bethe-Pierels approximation become exact on the Bethe lattice. 

6 Summary 

We have considered very general classical spin model with cubic symmetry on the Bethe lattice. 
First of all, for the underlying Cayley tree, which is the planar graph, we have presented the 
partition function of the model as a double power series in terms of graph expansions for 
arbitrary planar graphs. This double graph expansion representing the known formula for 
linear chain was evaluated exactly in case of Cayley tree. The expression for free energy per 
spin in the thermodynamic limit in this case was found to be continuous in temperature. This 
is in full agreement with general statement about mechanisms of phase transitions on recursive 
lattices. In the thermodynamic limit the influence of the huge amount of boundary sites, which 
in contrast to all other sites have only one neighbor, precludes the system from undergoing any 
phase transitions in the sense usually understood in the thermodynamics of " normal" systems. 
Then, getting rid of all boundary sites by formal putting of the coordination number of all sites 
equal to q, thus, passing to the completely homogeneous connected recursive graph, the Bethe 
lattice, we have applied the methods of the dynamical systems theory or, more precisely, the 
theory of discrete mappings. In this technique one exploits the self-similarity of the Bethe lattice 
and establishes a connection between the thermodynamic quantities defined for the lattices with 
different number of sites. We have used the method in which this connection is given by the 
RR for some formal variables x and y which have no direct physical meaning but can be used 
to determine any thermodynamic function of the system. The key point of this approach is 
the iterative procedure for the system of RR x n = f\ (x n _i, y n -i), Un — ji {x n -i, Un-i), where n 
can be regarded as the number of generations into the Bethe lattice. Thus, doing the iterations 
one approaches the thermodynamic limit n — > oo. However, actually the iterative procedure 
is terminated for some finite number of iterations by the fixed point (which is the case for 
ferromagnetic couplings)or by s-cycles, which takes place in the antiferromagnetic case. There 
is also more complicated case of the chaotic behavior of the iterative sequence but it is out of 
the topic of the article. We have identified the different thermodynamic phases of the system 
in ferromagnetic case (J > 0, K > 0) (disordered, partially ordered and completely ordered) 
with different types of the fixed points of RR. According to that we have obtained the phase 
diagrams of the model which are found to be different for Q < 2 and Q > 2. The former case 
contains three tricritical points whereas the latter just one tricritical and one triple points. 
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